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Abstract 

We calculate, to 3 loops in perturbation theory, the bare /3-function of QCD, 
formulated on the lattice with the clover fermionic action. The dependence 
of our result on the number of colors N , the number of fermionic flavors Nf, 
as well as the clover parameter csW) is shown explicitly. 

A direct outcome of our calculation is the two-loop relation between the 
bare coupling constant qq and the one renormalized in the MS scheme. 

Further, we can immediately derive the three-loop correction to the rela- 
tion between the lattice A-parameter and qq, which is important in checks of 
asymptotic scaling. For typical values of cswi this correction is found to be 
very pronounced. 
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I. INTRODUCTION 



The clover action for lattice fermions was introduced a number of years ago , as a means 
of reducing finite lattice spacing effects, leading to a faster approach to the continuum. It 
is widely used nowadays in Monte Carlo simulations of dynamical fermions. 

To monitor the onset of the continuum limit, tests of scaling must be performed on 
measured quantities. In particular, asymptotic scaling is governed by the bare /3-function, 
defined in the standard way: 



(1) 



(a is the lattice spacing, g (go) the renormalized (bare) coupling constant, /i the renormal- 
ization scale). In the asymptotic region go ^ one may write (3^ as: 



/3l{9o) = -bo 9o - bl gl - b^ gl + 
where, in SU {N) gauge theory with Nf fermion species, 
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(2) 

(3) 
(4) 



It is well known that the first two coefficients, bo and bi, are universal; the 3- loop coef- 
ficient 62 ) on the other hand, is regularization dependent. In the case at hand, 62 ^^us 
expected to depend not only on and Nf, but also on the free parameter csw which appears 
in the clover action (see next Section). 

In the present work we calculate 63 fo^' arbitrary N, Nf and csw- The analogous calcula- 
tion for the case of pure gauge theory without fermions, as well as for the case with Wilson 
fermions, was done a few years ago (Refs. and 0, respectively). We follow the general 
setup and notation of those publications, to which we refer the reader for further details. 

The /3-function enters directly into the relation defining the renormalization group in- 
variant parameter Al : 



aA] 



exp 



'2bog^ 
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l+qgl + 0(gt) 
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bl - bpb^ 
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(5) 



The "correction" factor q is, as we shall see, very pronounced for typical values of csw and 

A direct outcome of our calculation is the two-loop relation between the MS coupling 

a = g^/ (Air) and = 5'o/ (^tt) : 



a 



ao + di{ajj) al + d2{ajj) + O {o!^ 



(6) 



This relation is useful in studies involving running couplings or renormalized quark masses 
(see e.g. Refs. [^H^)- In fact, we will first compute this two-loop relation, which in turn will 
allow us to derive 62 ■ 
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In Sec. m we describe the method we employed to calculate 62- We present our results 
for di{afi) and d2{a^), as functions of N, Nf and csw- Discussion of further technical 
details and checks of our calculations is relegated to Sec. |ITT[ 

Given the difficulty of this computation, two different strategies have been adopted 
by each of the authors independently. H.P. performed a standard perturbation theory 
calculation, as described in Section II. A.B. used the approach of ref. [0, in which 
d2{afi) — d2{afi) |csw=o is obtained through a two-loop computation of the lattice Schrodinger 
functional. A detailed comparison of our results is presented at the end of Section III, show- 
ing very good agreeement. 



II. SETUP OF THE CALCULATION AND RESULTS. 

Our starting point is the Wilson formulation of the QCD action on the lattice, with the 
addition of the clover (SW) []T| fermion term. In standard notation: 

Sl = \ E Tr[l-f/^,,(x)]+EE(4^ + ^/.o)^/(a^)^/(a:) 

f X,fl 

1 _ ^ 

+ 7CSW E E i^f{x)a^^F^,^{x)ijf{x), (7) 

f X, ^, V 



where: F^,y = \[Q^y - Qy^^, Qp.v = U^,^y ^Uy^-y,^U-y,^-y ^U^^^^ (8) 

o 

Here Ufi^vix) is the usual product of link variables U^{x) along the perimeter of a plaquette 
in the /i-z/ directions, originating at x\ r is the Wilson parameter; / is a flavor index; m/,o 
are the bare fermionic masses; a^y = {i/2)['j^, ■j^]; powers of a have been omitted and may 
be directly reinserted by dimensional counting. The quantity cgw is a free parameter in the 
present work; it is normally tuned in a way as to minimize 0{a) effects. 

The procedure we have followed here is analogous to the one employed in Refs. @-§l; we 
shall only briefly describe it here, in the interest of a self-contained exposition. 

We set out to compute the relation between the bare lattice coupling and the renor- 
malized coupling g defined in the MS renormahzation scheme: 

9o = Zg{go, afi) g, (9) 

where /z indicates the renormahzation scale. Zg is related to the lattice /5-function (3^ {go), 
and to its MS-renormalized counterpart I3{g), through: 



f3{9) = ~3 l^'^^^^aido, ajj) \a,go= -bog^ - bi g^ - 626''^ + ••• (10) 
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The lattice /5-function, as defined above, is independent of tlie renormalized fermionic 
masses; tliis follows from the renormalizability of the theory |Q combined with dimensional 
arguments (see e.g. Ref. [Q). Thus, we will consider fermions with zero renormalized mass. 
From the above one can derive the relation 0] 

^\9o) = {l-gl^^\nZ^^ Z,i5{g^Z-^), (11) 

valid to all orders of perturbation theory. 

The one- and two-loop coefficients of Z| have the form: 

Zg(^o, a/i)^ = 1 + ^o(a/i) gl + Lx{a^) gl + 0(^^) ^^^^ 
Lo(a;) = 26olna; -|- Zq, -^i(a;) = 26i Inx -|- /i. 

The constant /q is related to the ratio of the A parameters associated with the particular 
lattice regularization and the MS renormalization scheme: 

/o = 26oln(AL/AMs)- (13) 

Its value as a function of N ^ Nf, csw; is known (see e.g. Ref. and references therein) and 
is presented here with increased accuracy for the csw-dependent coefficients: 

/o = -L _ 0.1699559991998031(2) N + Nf Iqi (14) 
where : loi{r=l) = 0.006696001(5) - cgw 0.00504671402(1) + cgw^ 0.02984346720(1) 

We have performed all our calculations for a range of r-values. In what follows we will 
be presenting only our results for r = 1, which is the value used in most Monte Carlo 
simulations; nevertheless, just for comparison, we give here /qi for r 7^ 1 : 

/oi(r=0.5) = 0.020173519(2) - csw 0.0144735322(1) + csw^ 0.0545958518(1) 
/oi(r=2.0) = -0.00181358(2) - cgw 0.0011966285(1) + cgw^ 0.0119267696(1) (15) 

We see that this range of r-values induces significant variations on Al; the dependence on 
Csw is also quite pronounced, as can be seen from Eq. (p!^, leading to changes in A^ of up 
to a factor of 2. 

Expanding Eq. (pH]) in powers of gQ leads to a tower of relations between the coefficients 
of P^{go) and P{g) ; namely, the difference between the (z-M)-loop coefficients, bf — bi , can 
be written in terms of i-loop contibutions to Zg. In particular, as is well known, b^ — bo = 0, 
bi — bi = 0, and, what will be important below: 

6^ = 62-61/0 + ^1- (16) 



Thus, since 62 is known from the continuum [|T0| 
1 



(An) 



2857,. / 1709A^2 i \ . fhQN 11 



54 M 54 36 4A^2 ) ^ / 27 18iV 



(17) 
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the evaluation of 62 requires only a two- loop calculation on the lattice, of the quantity h. 

We will compute Zg in the background field gauge |11-15|, where it is simply related to 
the background field renormalization constant Z^, 

ZAigo, a/i) Zg{go, afxf = 1. (18) 

In the background field formulation the links are written as 

where Qfi{x) = T^Q'^^^x) and = T'^A^^{x) are the quantum and background fields 

respectively; are Hermitian su{N) generators. We use the standard gauge-fixing term 
for this formulation, along with the corresponding Faddeev-Popov action for ghost fields, 
and the contribution from the integration measure (in our notation, these terms are written 
explicitly in [Q). 

As in Ref. 0], we write the renormalized one-particle irreducible two-point functions of 
the background and quantum fields as 



00 



1=1 

00 

r^'^(p)?. = -S'^' [iS,.p' - P.P.) (1 - ujM) + \p,p,] , u:n{p)= ^ g^'J^{p). (21) 

1=1 

Correspondingly on the lattice 

00 

i:^t\p)t = -^"'^P" [1 - ''(P)] /9l Kp)= T.9>^'\p)^ (22) 

M 1=1 

00 

E r?^(p)?, = -S'^'p' [3 (1 - uip)) + Ao] , u{p)= E g'o'uj'^^^p). (23) 
M 1=1 

There follows 

^^-^^-i-MpY 

The bare and renormalized gauge parameters appearing above, Aq and A, are related by 
A = ZqXq, where Zq is the renormalization constant of the quantum field: 

[l-u;n{p)] = ZQ[l-u;{p)]. (25) 

Both z/(p) and Lj{p) depend on (ap), go and Aq, whereas I'R^p) and ujr{p) depend on (p/fJ^), 
g and A. One can fix Aq = 1, and we have done so; this entails computing A = Zq, which is 
only needed here to one loop. 

The MS renormalized functions entering into Zg to two loops: i^r\p), i^R\p) and 
Vr {p) |a=i were calculated and presented in Ref. ||^ (and Ref. [|1^] for iVy = 0). 

The lattice two-point functions, u^^^p), uj'^^\p) and z/*^^^(p), were the main object of 
the present work, for csw 7^ 0. We will present here directly our results; details on our 
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calculations will be provided in the next section. For easier reference, we also include below 
results pertaining to pure gauge theory without fermions (from Refs. 0,^) and with fermions 
at csw = (from Ref. [g). 



z/(i)(p) 


Ao=l — 


z/(i)(p) 


Ao= 


--l,Nf=0 


+Nf 




Ao=l = 




Ao= 


=l,7V/=0 




z/(2)(p) 


Ao=l = 


Z/(2)(p) 


Ao = 


=1, Ar^=0 





^ ln(a V) + 



(26) 



1 1\, .00, , 1 



3N-- HaY) + ky- + hfN 



167r2 2 V NJ ' ^ ' N 



The Nf=0 quantities are 

UN . . o o. 1 



J^^'Hp) |Ao=i,7Vf=o = -— ^ln(aV) - ^ + iV 0.217098494366724... 

5iV . . o o^ 1 



cc^^^Hp) lAo=i,iv,=o =--—ln(ay)--^ + iV 0.137286278290915... (27) 

487r^ 8A' 



^/^'Hp) |Ao=i,iv,=o = -TjTT^ In(ay) + — — - 0.016544619540(4) + 0.0074438722(2). 
For r = 1 we found 



= - 0.013732194(5) + cgw 0.00504671402(1) - cgw 0.02984346720(1), 

k2f = 0.0011877(14) - Csw 0.0001578(3) - csw^ 0.0052931(2) 

- csw^ 0.00050624(3) - csw^ 0.00008199(1), (28) 

k^j = - 0.0013617(16) - Csw 0.0000981(4) + csw^ 0.0052440(6) 

+ csw^ 0.00021431(3) + csw^ 0.00004382(1). 

The origin and meaning of the errors in the above results will be explained in the next 
section. Contributions which vanish in the continuum limit, a — > 0, do not contribute to Zg 
and are neglected as in Refs. pj-Q. 

We have now what is needed to calculate Zg to two loops using Eq. (^if) . The resulting 
two-loop constant li is given by 



h = ^ + 0.018127763034(4) - 0.0079101185(2) + 



128A^2 



/ 



^11 — + /l2 

N 



l^^ = - 0.0011967(14) + csw 0.0001578(3) + csw^ 0.0052931(2) 

+ csw^ 0.00050624(3) + csw^ 0.00008199(1), 

= 0.0009998(16) + csw 0.0000342(4) - csw^ 0.0048660(6) 

- csw^ 0.00021431(3) - csw^ 0.00004382(1). 



(29) 



(30) 



At this point, it is a trivial exercise to recover 63 ) fo'^ ^iny value of A^, Nf, csw, by 
substitution in Eq. (0). We list here some particular cases of interest, setting = 3: 

6i'(csw=0) = - 0.00159983232(13) + Nf 0.0000799(4) - N] 0.00000605(2) 

h^[cs^=l) = - 0.00159983232(13) - Nf 0.0009449(4) + Nj 0.00006251(2) . 

b^(Nf=3) = - 0.0014144(10) + csw 0.0000654(2) - csw^ 0.0024241(4) 

- csw^ 0.00008108(2) - csw^ 0.00001781(1). 
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The correction coefficient q of Eq. also follows immediately. In the particular cases 
considered above we obtain (A^ = 3): 

csw=0 : q = 0.18960350(1) (iV/=0), q = 0.2160(1) (A^/=2), q = 0.2355(2) (A^/=3) 

csw=l : q = 0.18960350(1) (iVj=0), q = 0.4529(1) (A^/=2), q = 0.6138(2) (A^/=3) 

q{Nf=3) = 0.2355(2) - cgw 0.01007(4) + cgw^ 0.3731(1) + cgw^ 0.01248(1) + cgw^ 0.002741(1). 

(32) 

Clearly, these values of q bring about a substantial correction to asymptotic scaling, with a 
pronounced csw dependence. 

Finally, the coefficients di^a/i) and d2{an), in the relation (^ between the MS coupling 
a and a^, can be read off Lo{x) and Li{x), defined in Eq. (p!2D: 



di{x) = -47rLo(x), d2{x) = (Airf Loixf - Li{x) (33) 



III. THE CALCULATION IN LATTICE PERTURBATION THEORY 

We describe here briefly some technical features and checks of our calculation of the 
quantities u^^\p), u^'^^p). Fermionic contributions to iu^^\p) are identical to those in u^^^p). 
Many aspects are as in Ref. |Q, and are discussed there in more detail. 

Two diagrams containing fermions contribute to iy^^\p), shown in Figure 1. There are 
18 two- loop fermionic diagrams contributing to u^'^\p), as well as two diagrams containing 
an insertion of the one-loop fermion mass counterterm; these are shown in Figure 2. 




1 2 

FIG. 1. Fermion contributions to the one-loop function i'^^\p). Dashed lines ending on a cross 
represent background gluons. Solid lines represent fermions. 

The algebra involving lattice quantities was performed using a symbolic manipulation 
package in Mathematica, developed by one of the authors (H.P.) and collaborators over the 
past years (see @] and references therein). In order to carry out the present work, this 
package was augmented to include vertices of the clover action. These vertices are more 
complicated than those of Wilson fermions, as can be seen by the form of the action (|^); as 
a result, the algebraic expression for each diagram is considerably longer in the present case. 
On the other hand, all new contributions to fermionic vertices are accompanied by extra 
powers of momenta, and consequently lead to expressions with improved infrared behaviour. 

The two-loop amplitude u^"^^ (p) can be written as: 
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.(2) 



(34) 



(the index i runs over the diagrams shown in Fig. 2). Since two-loop diagrams can have a 
maximum of 4 fermionic vertices, t'i(p) is a polynomial in csw of degree up to 4, starting 
from csw^- 

On general grounds we expect: 



(p2 — 4^^sin^(p^/2)), where the most general dependence of Cn,i on N, Nf, csw is: 



+ CM^ + C44+0((apr) (35) 



AT, 



(36) 




5 6 7 8 




9 10 11 12 




13 14 15 16 




17 18 19 20 

FIG. 2. Fermion contributions to the two-loop function u^'^\p). Dashed lines represent gluonic 
fields; those ending on a cross stand for background gluons. Solid lines represent fermions. The 
filled circle is a one-loop fermion mass counterterm. 
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As a byproduct of improved infrared convergence, the coefficients ci^j, C2,i vanish for each 
diagram separately in this case; all other coefficients are also separately infrared convergent. 

To extract the p-dependence, we apply a procedure of successive subtractions ^ on terms 
containing (sub-) divergences; our expressions then split into a relatively small part made of 
factorized one-loop terms, whose (possibly logarithmic) p-dependence is easily extracted in 
terms of tabulated lattice integrals, and a much larger part which can be evaluated by naive 
Taylor expansion in ap. 

Once the p-dependence has been made explicit, the resulting expressions for c„ » must 
be numerically integrated over internal momenta. For each expression, our package creates 
highly optimized numerical code for integration in momentum space, over finite lattices of 
size L . An extrapolation to infinite size is then performed. This operation is of course a 
source of systematic error. In a nutshell, we estimate this error as follows: First, different 
extrapolations are performed using a broad spectrum of functional forms of the type: 

Ee„L-(lnL)^- (37) 

For the k^^ such extrapolation, a deviation is calculated applying three different quan- 
titative criteria for quality of fit. Finally, these deviations are used to assign weights 
dk'^ / {Y,k dk'^) to each extrapolation, which are then combined to produce a final value to- 
gether with the error estimate. We have checked the validity of these estimates in cases 
where the exact answer was known, finding the estimate to be always correct. 

The improved infrared behaviour of each diagram led to a rather stable extrapolation; 
indeed, for the accuracy reported in this work, lattices up to L < 28 were sufficient in all 
but a few cases. 

The results for each diagram are presented in Tables I through VI. Diagrams not ap- 
pearing in these Tables give vanishing contributions. 

For the one- loop amplitude op^ u^^\p), the individual contributions from the correspond- 
ing two diagrams in Figure 1 are: 

u^(p) = csw 0.00504671402(1) - cgw^ 0.02984346720(1), z/2(p) = (38) 

There exist several constraints on the coefficients Cn,j, which we have used as verification 
both on the algebraic expressions and on the numerical results: 

> Co,i: Gauge invariance requires J2iCo,i = 0. This relation amounts to 4 independent 
nontrivial checks, for the coefficients of csw csw^ Csw/A^; Csw^/A^- Our algebraic 
expressions pass these checks. Our numerical results also give zero, well within the 
systematic errors; this effectively controls the consistency of our error estimation. 

> C3_j: The total contribution for single logarithms must give zero: = 0. Again, 
this relation amounts to 5 independent nontrivial checks, which were performed both 
algebraically and numerically. 

> There are several other relations among the coefficients appearing in individual dia- 
grams (e.g. C3^io^ = 4c3^'^3'', c'^'iy = —5c^'ir^, etc.); these can be understood in terms 
of the vertices entering these diagrams, and have been used as further consistency 
controls. 
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The most stringent check is provided by Ref. in which the Schrodinger functional in 
lattice QCD is computed to two loops, using the clover action with csw = 0, 1. Eq. (5.6) of 
that reference^] reads: d2i{N=3, Csw=l) = 1-685(9) — 8.6286(2) Cg^. Calculating the same 
quantity from our results, we deduce: 1.6828(8) — 8.62843775(1) Cg\^- Both sets of numbers 
are clearly in very good agreement. 

Since the dependence of d2 on cgw is polynomial of degree 4, all coefficients of the 
polynomial (except for the 0th degree, which is gotten from |^) can be checked by repeating 
the calculations of Ref. for a few more values of csw- We have done so for all values of cgw 
listed in Table VII; the results obtained in this way have a larger error than those stemming 
from Eqs. (|33| , |T2|, |2^, jS^), but they provide a very solid and thorough confirmation of our 
main result, Eq. (PPD . 

The entire procedure employed in this paper can be applied without modifications to a 
number of other cases of interest, among them the higher order renormalization of 0{a)- 
improved operators. We hope to study these cases in a future work. 

Acknowledgements: H.P. would like to thank P. Weisz for many helpful discussions and 
suggestions. 



^In the notation of Ref. 0: csw = + Cg^ + ■ • • 



sw ^ ^sw ■ 
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TABLE I. Coefficients 4 r = 1. 



J ^;Oj_i 

1 0.000040200(3) 

12 0.00114029(4) 0.00114886(4) 

14 -0.00008041(1) 

15 -0.00013096(2) -0.00007625(1) 

18 -0.0009691(2) -0.0010726(1) 

19 0.002862565(1) 0.001191276(1) 

20 -0.0028625645(1) -0.0011912762(1) 



TABLE IL Coefficients c^r^'/K r = 1. 



i 




Jl>2) 


4 


-0.00016080(1) 


0.0071556(1) 


5 


0.000040202(2) 


-0.0077678(3) 


11 


-0.00021503(4) 


0.0023875(1) 


12 


-0.00114029(4) 


-0.00114886(4) 


13 


-0.0000575(1) 


0.000956(1) 


14 


0.00029874(1) 


0.00026811(2) 


16 


0.00004358(1) 


-0.00011562(2) 


17 


0.0002220(2) 


-0.002807(1) 


18 


0.0009691(2) 


0.0010726(1) 


19 


-0.002862565(1) 


-0.001191276(1) 


20 


0.0028625645(1) 


0.0011912762(1) 



TABLE III. Coefficients 4 l'^^- r = 1. 



i c^~^'^^ c^".^'^^ 

"8 -0.30986678046212 

14 -0.01449434463(1) 

15 0.00500026042(1) -0.00589886486(1) 

18 -0.0933193205(1) 0.0566733951(1) 0.03619153756(1) 

19 0.41268018519(1) -0.05077453020(1) -0.03619153756(1) 
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TABLE IV. 



Coefficients c^'P- r 



10 
13 

14 
16 
17 
18 
19 



.(1.1) 

-3,i 



0.30986678046212 

0.01009342804(1) 

0.00252335701(1) 

-0.00500026042(1) 
0.0271111297(1) 

-0.02523357012(1) 
0.0933193205(1) 

-0.41268018519(1) 



„(1.2) 





-0.11937386878(1) 
-0.02984346720(1) 

0.01052896544(1) 
-0.00463010057(1) 

0.14921733598(1) 
-0.0566733951(1) 

0.05077453020(1) 



J1.3) 










-0.03619153756(1) 
0.03619153756(1) 



TABLE V. Coefficients c[J'^\ r = 1. 



i 


J-i.i) 


J-1.2) 


J-1.3) 


J-1.4) 


3 


0.00082631538(1) 


-0.007460866799(1) 








4 


0.000020100(1) 


-0.00119261(2) 
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0.003071859(3) 


0.00121700384(1) 








11 


-0.000020100(1) 


0.00059631(2) 








14 


0.00018663(2) 


0.000005901(1) 








15 


0.0000931(2) 


0.0001160(1) 


-0.00000655(2) 


-0.000038170(1) 


18 


0.0004905(4) 


-0.0002038(3) 


-0.00059685(2) 


0.00008405(1) 


19 


-0.004826287(1) 


0.0016290054(5) 


0.0000971639(1) 


-0.0001278777(1) 



TABLE VI. Coefficients c^^f. r = 1. 



i 


^(1,1) 


J1.2) 


J1.3) 


Jl>4) 


3 


-0.00093961433(1) 


0.008800841023(1) 








4 


-0.0000050250(3) 


0.000298152(4) 








5 


0.0000004358(3) 


-0.0001706(1) 








8 


-0.003071859(3) 


-0.00121700384(1) 








10 


-0.000169546987(1) 


0.002005213655(1) 








11 


0.00003775(1) 


-0.00074560(4) 








13 


0.00004011(2) 


-0.0005579(2) 








14 


-0.0002223(1) 


-0.0001602(1) 


-0.00001405(1) 





16 


-0.0001813(4) 


-0.0000980(5) 


-0.00027132(2) 





17 


0.0000775(3) 


-0.001486(1) 








18 


-0.0004905(4) 


0.0002038(3) 


0.00059685(2) 


-0.00008405(1) 


19 


0.004826287(1) 


-0.0016290054(5) 


-0.0000971639(1) 


0.0001278777(1) 
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TABLE VII. Comparison of our results for {d2i 
dard perturbation theory (PT) (Eqs. (||, |l|, |2|, | 
functional (SF) methods of Ref. [7]. {N = 3, 4\J, 



-d2i I (0) _ )/(167r^) coming from stan- 
to those we obtain using the Schrodinger 



0) 



.(0) 



-4.0 
-2.0 
-1.0 
-0.5 
0.5 
1.0 
2.0 
4.0 



PT 



0.202270(29) 

0.049518(10) 

0.012619(7) 

0.003233(2) 

0.003197(2) 

0.013257(7) 

0.056484(10) 

0.261725(29) 



SF 



0.2029(21) 

0.04968(53) 

0.013041(60) 

0.003258(96) 

0.003191(18) 

0.01327(4) 

0.05659(45) 

0.2627(33) 
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